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Biermann Mechanism in Primordial Supernova Remnant and Seed 
Magnetic Fields 
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Abstract. We have studied the generation of magnetic fields by the Biermann mechanism in the pair-instability supernovae 
explosions of the first stars. The Biermann mechanism produces magnetic fields in the shocked region between the bubble 
and interstellar medium (ISM), even if magnetic fields are absent initially. We have performed a series of two-dimensional 
magnetohydrodynamic simulations with the Biermann term and estimate the amplitude and total energy of the produced 
magnetic fields. We find that magnetic fields with amplitude lO^^** — 10~^^ G are generated inside the bubble, though the 
amount of magnetic fields generated depend on specific values of initial conditions. This corresponds to magnetic fields with 
total energy of 10^* — 10"'^ erg per each supernova remnant, which is strong enough to be the seed magnetic field for a 
galactic and/or interstellar dynamo. 
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1. Introduction 

Magnetic fields are ubiquitous in the universe. In fact, obser- 
vations of rotation measure and synchrotron radiation have 
revealed that magnetic fields exist in astronomical objects 
with various scales: galaxies, clusters of galaxies, extra- 
cluster fields, etc. (for a review, see, e.g., Widrow 2002). 
The observed galactic magnetic fields have both coherent and 
fluctuating components whose strengths are comparable to 
each other (Fosalba et al. 2002; Han et al. 2004). Conven- 
tionally, these magnetic fields are considered to be amplified 
and maintained by a dynamo mechanism. The coherent com- 
ponent in a galaxy is expected to be amplified by a galac- 
tic dynamo, while the fluctuating component may be ampli- 
fied by an interstellar dynamo driven by turbulent motion of 
the interstellar medium (ISM; Balsara et al.2004). However, 
the dynamo mechanism itself cannot explain the origin of the 
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magnetic fields: seed magnetic fields are needed for the dy- 
namo mechanism to work. 

Thus far, various mechanisms have been suggested as a 
possible source of the seed field. They can be classified into 
two types, those of astrophysical origin and those of cosmo- 
logical origin. Here we concentrate on the former (for the 
latter scenario, see, e.g., Davies & Widrow 2000; Takahashi 
et al. 2005; Amjad & Robert 2005; Tsagas 2005; Yamazaki 
et al. 2005). Basically, the astrophysical magnetogenesis in- 
vokes the Biermann mechanism (Biermann 1950), which is 
induced by the electric currents produced when the spatial 
gradient of the electron pressure is not parallel to that of 
the density. This is a pure plasma effect so that there is no 
need to assume unknown physics as is often done in cosmo- 
logical models. Because the Biermann mechanism requires 
a non-parallel spatial gradient of the pressure and density, 
some nonadiabatic process is necessary to produce devia- 
tion from a polytropic equation of state. The strong magnetic 
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fields in high-redshift galaxies (Athreya et al. 1998) imply 
that a significant amount of the seed magnetic field should 
be generated at an early stage, e.g., the epoch of cosmo- 
logical reionization or protogalaxy formation. For instance, 
Gnedin et al.(2000) studied the generation of magnetic fields 
in the ionizing front and found that magnetic fields as high as 
w 10~^^ G in virialized objects can be generated. Kulsrud et 
al.(1997) showed that a magnetic field of « 10^^^ G can also 
be generated at shocks of large-scale structure formation. 

In this paper, we investigate magnetogenesis at smaller 
scales. Specifically, we study the amplitude of the mag- 
netic field produced by the Biermann mechanism when the 
shock waves of the supernova explosions of the first stars are 
spreading throughout the ISM (Hanayama et al. 2005). The 
primordial supernova explosions are expected to take place 
effectively, since the initial mass function (IMF) of Popula- 
tion III stars should be substantially top-heavy (e.g., Abel et 
al. 2002). We perform a series of two-dimensional magne- 
tohydrodynamic (MHD) numerical simulations in which the 
Biermann term is included. We also discuss whether they can 
be the origin of the cosmic magnetic fields. Consequently, 
we find that the spatially averaged amplitude of the produced 
magnetic field in virialized objects reaches 10^^^ G, which 
is much greater than those expected from cosmic reioniza- 
tion and large-scale structure formation. Thus, the supernova 
explosions of the first stars can be effective sources for the 
seed magnetic fields. Although the situation considered here 
is somewhat similar to that of Miranda et al.(1998), they as- 
sumed a multiple-explosion scenario of structure formation 
and considered explosions of objects with mass > IO^Mq 
at z > 100, which is clearly unrealistic in the context of the 
current standard model of structure formation. 

2. Numerical simulations and results 

We solve the two-dimensional MHD equations using con- 
served quantities with heating and cooling adjusted for the 
ISM. The induction equation with Biermann term is as fol- 
lows: 

Vp X VP 



- = Vx(vxB) 



(1) 



where p, v, P, and B are the density, velocity, pressure, and 
magnetic field, respectively. The last term of the right-hand 
side of equation Q is the Biermann term: a in equation 
Q is the so-called Biermann coupling constant defined by 
a — TOpc/e(l + x) ^ 10^* G s, where nip, e, and x are 
the proton mass, electric charge, and ionization fraction, re- 
spectively. Although the gas temperature just before star for- 
mation begins is rather cool (T 200K; Abel et al. 2002; 
Bromm et al. 2003; Omukai & Palla 2003), UV radiation 
from the first stars ionizes the surrounding ISM (Freyer et 
al. 2003; Mori et al. 2004). Thus we set x — 5/6 assuming 
n-He/nn '^0.1. We use a cooling function derived by Ray- 
mond et al.(1976). When parcels cool below 10'* K, an artifi- 
cial heating rate proportional to the density is used. The con- 
stant heating coefficient is set so that heating balances cooling 
at the ambient density and temperature. Although the cooling 
function and heating rate in the primordial gas are not clear so 



far, they are not important for the adiabatic expansion phase 
that we concentrate on. 

We solve the above equations by the two-dimensional 
MHD code in cylindrical coordinates (r, z, cj)) assuming axial 
symmetry around the symmetry axis (z). The code is based 
on the modified Lax-Wendroff scheme with an artificial vis- 
cosity of von Neumann and Richtmyer to capture shocks. The 
numerical scheme was tested by comparing known solutions 
that have been obtained either analytically and numerically 
(see e.g., Hanayama & Tomisaka 2005). 

In all the computations, grid spacings are chosen Ar = 
Az = 0.1 pc. For example, the numerical domain covers a 
region of 130 pc x 130 pc with 1300 (r) x 1300 (z) mesh 
points in our fiducial model. We begin the simulation by 
adding a thermal energy of Eq = 10^^ or 10^^ erg within 
the sphere of 2 pc in radius. 

As the bubble expands, the ejected gas interacts with the 
ISM, and a shock wave is formed. In the shocked region the 
gas is heated nonadiabatically, which is a necessary condi- 
tion for the Biermann mechanism to work. Here the structure 
of the interstellar environment is important because it affects 
the density and pressure profiles of the shocked region that is 
directly related to the Biermann term. We assume an inhomo- 
geneous ISM with average density tiism = 0.2 cm^'^. This 
is roughly consistent with the situation discussed in Bromm 
et al.(2003). The scale length of the density and the ampli- 
tude of the inhomogeneity are as yet poorly understood, and 
we assume inhomogeneity with the scale length A = 1 pc 
and density variation 0.2 x 2*^cm^'^, which are values sim- 
ilar to those in our galaxy. Within the variation, the ampli- 
tude of the density is given at random and the distribution 
is smoothed numerically to create a perturbation. This is our 
fiducial model for the ISM. As we show in the next section, 
the amplitude of the produced magnetic field is sensitive to 
the scale length A, while the average density and the ampli- 
tude of the density variation are rather unimportant. Thus we 
consider several different models in addition to the fiducial 
model: specifically, we vary the mean density [1 x 2*^ cm~^ 
and 10 X 2*^ cm^-^] and the scale length (3 and 10 pc). 

As for the explosion energy of the supernova, we adopt 
^-SN = 10^'^ erg for the fiducial mode. This explosion en- 
ergy corresponds to stars with mass 250AfQ that explode as a 
pair-instability supernova (Fryer et al. 2001). In addition, we 
consider a model with Esn = 10^^ erg as a variation. 

Figure ^ shows the contours of the gas density and the 
amplitude of the magnetic field produced by the end of the 
adiabatic expansion phase t = 1.26 x 10^ yr for the fiducial 
model. The radius of the bubble is about 125 pc, and turbulent 
motion is induced in the shocked region because of the inho- 
mogeneity of the ISM. The amplitude of the magnetic field 
is about 10~^* G for the central region and about 10^*^ G 
just behind the shock. The total magnetic energy inside the 
bubble is about 10"^° erg. 

We have also checked the case of a homogeneous medium 
to test the robustness of our computation. We find that, on 
average, a SNR generates magnetic fields with an amplitude 
^ 10^^^ G behind the shock front, which is ^ 10^"^ times 
smaller than that in an inhomogeneous medium. Therefore, in 
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the case of a homogeneous medium, we estimate a numerical 
error of the amplitude of the magnetic field of 10^^^ G. It 
should be noted that the resulting magnetic field has only a 
toroidal component because axial symmetry was assumed. 

In Figure|2] we show the time evolution of the total mag- 
netic energy for various models. The behaviors are qualita- 
tively similar for all the models. The robust knees around 10'^ 
yr in Figure|2]come from the formation of an adiabatic shock 
front: this corresponds to when the SNR shifts from the free 
expansion phase to the Sedov phase. The total magnetic en- 
ergy at the end of the adiabatic expansion phase is larger for 
models with a smaller scale length of the ISM density. This 
tendency is confirmed by an order-of-magnitude estimate in 
the next section. For models with large average ISM den- 
sity or small explosion energy, the total magnetic energy is 
smaller because the bubble is smaller than in the other mod- 
els. Although there are many uncertainties in the initial con- 
ditions, the generation of a magnetic field with a total energy 
of 10^* — 10"^^ erg appears to be robust. 

3. Analytic estimates of magnetic fields and 
implications for the seed magnetic Field 

To understand the result of the numerical simulations, in this 
section we perform an order-of-magnitude estimation of the 
strength of the magnetic field produced by the Biermann 
mechanism. The amplitude of the magnetic field produced 
by the Biermann mechanism can be estimated from the Bier- 
mann term in equation Q, 

Vp X VP 
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where At ~ 10'^ yr is the characteristic timescale in which 
the Biermann mechanism works. Taking the characteristic 
pressure to be the ram pressure of the gas, P ~ Pram = 
p f^ub' '^^'^ '^he characteristic velocity of the bubble vt,uh ^ 
10~^ pcyr^^, we obtain 
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where L is the scale length of the pressure component per- 
pendicular to the density gradient. Then the magnetic energy 
produced by the Biermann mechanism for each primordial 
SNR can be estimated as 
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which is consistent with the value obtained from our numeri- 
cal simulations. 



The dependence of the total magnetic energy on several 
parameters can also be understood from equation ©. It is 
found directly from equation that Eb oc A^^. To examine 
the dependence on the other parameters, we simply assume 
the Sedov-Taylor solution: 

1/5 / P \ 1/5 

,"-ISM/ 



Pbub tX t 



2/5 



-E'SN \ ,-3/5 



(5) 



Putting these into equation yields Eb oc (-E'sn/'t-ism)^^^- 
Our numerical results shown in Figure|2]are quite consistent 
with these simple estimations. 

Now we estimate the spatially averaged energy den- 
sity of the magnetic fields produced by the first stars and 
consider whether they can be a source of the seed fields. 
For the primordial star formation rate, we extrapolate the 
one by Pello et al.(2004) and Ricotti et al.(2004); p^, ^ 
10~^ Mq yr~iMpc^'^. Denoting the magnetic energy pro- 
duced by the Biermann mechanism as €sn ^ lO'^^erg, the 
magnetic energy density produced during the formation pe- 
riod of the first-star (r ^ 1 Gyr) can be obtained as 
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where Msn is the typical mass scale of first stars that end 
up in pair-instability supernovae, and /^^ is the mass frac- 
tion of such first stars; we adopt j-^^ — 0.06, which was de- 
rived under the assumption that very massive black holes pro- 
duced from the first stars end up in supermassive black holes 
in galactic centers (Schneider et al. 2002). We note that the 
value is the comoving density averaged in the universe; we 
can convert the value to physical density in virialized objects 
(i.e., protogalaxies) as 
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where A is the density contrast. This corresponds to the 
mean magnetic field of i? IQ-^^ G in protogalaxies, 
which is much stronger than that expected in ionizing fronts, 
B - 10-1^ G (Gnedin et al. 2000). 

4. Summary and discussion 

We have studied the generation of magnetic fields in 
primordial supernova remnants. We have performed two- 
dimensional MHD simulations with the Biermann term, 
which can produce a magnetic field through the nonadiabatic 
interaction between the bubble and ISM, even if there is no 
magnetic field at first. We have found that the ISM around 
the primordial supernovae is an effective site for producing 
magnetic fields. The total energy of the magnetic fields is 
Xq28 _ ]^q31 gj,g^ depending on the parameters adopted. On 
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Fig. 1. Contours of the gas density (left) and the magnetic 
field (right) at the end of the adiabatic expansion phase 
(t = 1.26 X lO'^yr) for the fiducial model. The amplitude 
of the magnetic field is about 10^^* G at the central region 
and about 10^^^ G just behind the shock. 



the basis of the results, we have estimated the spatially av- 
eraged energy density of the magnetic fields produced by the 
first stars during the formation period of the first stars. The av- 
eraged energy density is about 10^'*" erg cm~^, which cor- 
responds to i? ~ 10^^^ G in protogalaxies at z ~ 10. This 
is much greater than expected from cosmic reionization and 
large-scale structure formation. Thus primordial supernova 
remnants would be a promising source for the seed fields for 
the galactic and/or interstellar dynamo. 

Although the coherence length of the seed field computed 
here is much smaller than the galactic scale, it can be am- 
plified by the galactic dynamo to produce a coherent com- 
ponent if the coherence length is about 100 pc (Poezd et al. 
1993; Ferriere & Schmitt 2000), which is a typical size of 
supernova remnants. It might also be amplified by an inter- 
stellar dynamo to produce the fluctuating component (Bal- 
sara et al. 2004). While it is beyond our scope to discuss the 
relation between the magnetic field produced and the galac- 
tic/interstellar dynamo, we plan to investigate the evolutions 
of the seed magnetic fields computed here as a result of the 
dynamo processes on the large scale of galaxies and clusters 
of galaxies. This will be presented in a forthcoming paper (H. 
Hanayama et al. 2005, in preparation). 

Regarding their use as an observational signature, a pro- 
posal to detect seed magnetic fields was made by, e.g., 
Plaga(1995). If there exist intergalactic magnetic fields pro- 
duced by the primordial SNRs, the arrival time of high energy 
gamma-ray photons from extragalactic sources would be de- 
layed by the action of intergalactic magnetic fields on elec- 
tron cascades (Ando 2004). Even a magnetic field as weak 
as ~ 10^^"' G would be detectable if the delay of the arrival 
time comes within a reasonable range (a few days). Therefore 
the mechanism of magnetic field generation proposed in this 
paper might be tested by the future high-energy gamma-ray 
experiments such as GLAST (Gamma-Ray Large Area Space 
Telescope). 
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Fig. 2. Time evolution of the total magnetic energy produced 
by the Biermann mechanism for various models. 
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